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In a viscous lock-exchange gravity current, which describes the reciprocal exchange of 
two fluids of different densities in a horizontal channel, the front between two Newtonian 
fluids spreads as the square root of time. The resulting diffusion coefhcient reflects the 
competition between the buoyancy driving effect and the viscous damping, and depends 
on the geometry of the channel. This lock-exchange diffusion coefficient has already 
been computed for a porous medium, a 2D Stokes flow between two parallel horizontal 
boundaries separated by a vertical height, H, and, recently, for a cylindrical tube. In the 
present paper, we calculate it, analytically, for a rectangular channel (horizontal thickness 
b, vertical height H) of any aspect ratio {H /h) and compare our results with experiments 
in horizontal rectangular channels for a wide range of aspect ratios (I/IO — 10). We also 
discuss the 2D Stokes-Darcy model for flows in Hele-Shaw cells and show that it leads 
to a rather good approximation, when an appropriate Brinkman correction is used. 



1. Introduction 

The lock-exchange configuration refers to the release, under gravity, of the interface 
between two fluids of different densities, confined in the section of a horizontal channel. 
This physical process has prompted renewed interest, as a part of the carbon dioxide 
sequestration issues (Neufeld & IIuppert(2009) ). The top of Fig. [I] shows the initial lock- 



exchange situation of a so-called full-depth release. The two fluids, initially separated 
by a vertical barrier (the lock gate), flll the whole section of the tank. When the gate 
is withdrawn (bottom of Fig. [T]), buoyancy drives the denser fluid along the bottom 
wall, while the lighter one flows in the opposite direction at the top of the channel. 
The so-called lock-exchange results in the elongation of the interface between the two 
fluids along the horizontal direction. Different regimes have been reported for the ve- 
locity and shape of the elongating interface. The slumping phase refers to the initial 
regime where inertia dominates over viscous forces, which typically applies for the case 



of salted and fresh water in a tank. In this regime, Benjamin(I968) , and more recently 



Shin et aZ.(2004)Shin, Dalziel & Linden showed that in the presence of a small density 
contrast (i.e. in the Boussinesq approximation Ap << p), the two opposite currents 
traveled at the same constant velocity. When the Boussinesq approximation does not 



apply, Lowe et a/.(2005)Lowe, Rottman & Linden 


Birman et aL(2005)Birman, Martin 


& MeiburgI Cantero et a/.(2007)Cantero, Lee, Balachandar & Garcia and 


Bonometti 


et aZ.(2008)Bonometti, Balachandar & Magnaudet 


showed that the two opposite fronts 



did travel at constant, but with different velocities. However this interface elongation, 
proportional to the time, is slowed down at later stages, in the viscous phase, where 
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dissipation prevails over inertia. In the latter regime, the interface elongates as t", where 
the exponent a, smaller than unity, may take different values depending on the geometry 
and the confinem ent of the flow ([Pidden fc Maxworthy(1982) Huppert(1982)[ Gratton 



& Minotti(I990) ; 'Cantero et a/.(2007)Cantero, Lee, Balachandar fc Garcia]^ Takagi fc 



Huppert(2007) Hallez fc Magnaudet(2009)[ ). In porous media, |Bear(1988)| and [Huppert, 
fc Woods(1995)| predicted an interface spreading proportionally to t he square root of 



time that |Seon et aL(2007)Seon, Znaien, Salin, Hulin, Hinch fc Perrm observed in a hor- 
izontal cylindrical tube. Such a spreading can be quantified with a diffusion coefficient, 
which reflects the balance between the buoyancy driving and the viscous damping. This 
coefficient, which depends on the nature and the geometry of the flow, has been com- 
puted for a porous medium by Huppert fc Woods(f995) for a 2D Stokes flow between 
two parallel horizontal boundaries separated by a vertical height, H, by |Hinch(2007) and 



Taghavi et aL(2009)Taghavi, Seon, Martinez fc Frigaard[ and for a cylindrical tube by 



Seon et aL (2007) Seon, Znaien, Salin, Hulin, Hinch fc Perrin However, to our knowledge, 



such a diffusion coefficient has not been derived for a rectangular channel (horizon- 
tal thickness 6, vertical height H, Fig. [T]), for which one expects to recover the porous 
medium regime for 6 <C -ff , and, possibly, the Stokes flow regime for 6 ^ iJ. In order to 
gather the limiting cases in the same paper, we first recall the results for porous media 
and 2D Stokes flows, together with the tube case, for the sake of comparison. Then we 
compute, for a rectangular channel of aspect ratio, the dependence of the interface 
h{x, t) and the corresponding viscous lock-exchange diffusion coefficient. We also test the 
so-called Stokes-Darcy 2D model to this lock-exchange configuration. Finally, we test 
and validate our theoretical results with experiments in horizontal rectangular channels 
for a wide range of aspect ratios (f/fO — 10). 



2. Lock-exchange in different geometries 

Let us first recall the basic hypotheses on the viscous gravity currents, common to 



the different geometries, used for instance by Huppert & Woods(1995) or Hinch(2007) |. 
As sketched in Fig. [l] the interface between the two fluids is assumed independent on 
the y direction. Its distance from the bottom boundary of the vessel is denoted h{x,t). 
This interface can be a pseudo-interface between two miscible fluids for which molecular 
diffusion can be neglected or between two immiscible fluids, provided that the interfacial 
tension can be neglected. The flow is assumed to be quasi-parallel to the horizontal x 
axis. This is a key hypothesis. Neglecting accordingly the vertical component of the fluid 
velocity implies that the vertical pressure gradient follows the hydrostatic variations: 
dp/dz = —pg. This hypothesis is violated at short times, immediately after the opening 
of the gate, but should become valid at later stages, as soon as the interface has slumped 
over a distance larger than H, thus ensuring a small enough local slope dh/dx. Then, 
the pressures, P+ and P_, in the lower layer, < z < h(x,t), and in the upper one, 
h{x,t) < z < H, respectively, write 

P+ = p{x, t) - p+gz ; P_=P+ + Apg{z - h{x, t)) (2.1) 

where p{x, t) denotes the pressure at the lower wall, z = 0. The difference between the 
horizontal pressure gradients in the two fluids is therefore linked to the interface slope 
by: 

dP+ aP_ , dh 
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Figure 1 . Sketch of the rectangular ceU of height H in the gravity direction and width b. Top: 
Initial configuration where the heavy fluid of density, p, is separated from the light fluid of 
density, p — Ap, by a vertical gate (plane y — z). Bottom: After removing the gate, buoyancy 
differences cause the denser fluid to flow in one direction along the bottom of vessel, while the 
lighter one flows in the opposite direction at the top of the vessel. The goal is to determine the 
space and time dependencies of the pseudo-interface, h{x,t). 



The time evolution of the interface, h{x,t), is governed by the mass conservation of each 
fluid (see Fig. [l] for notations). For instance, for the heavier bottom layer we have: 

dh dq 

where q(x,t) — q^{x,t) is the horizontal flux {ni? / s) of the denser fluid at the location 
x: 



q{x,t) = q+{x,t) = 



1 



u^{x,y,z,t)dydz 



(2.4) 



with Ux{x, y, z, t) the x— velocity component and b the spanwise length (the y integration 
is along this spanwise length). Moreover, in our configuration of uniform section along 
the horizontal axis x, Ux{x,y, z,t) must also satisfy the no net flux condition: 



q+ix,t) +q-ix,t) 



H 



Ux{x,y, z,t)dydz = 



(2.5) 



We will see in the following that, in the viscous regime of interest, the horizontal velocity 
component Ux, solution of either a Darcy or a Stokes equation, is proportional to the 
pressure gradient in each fluid layer. Such solutions, combined with eq. (2.2), eq. (2.4 1 
and eq. (2.5) allow then to eliminate the pressure gradients and to derive an expression 
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of the flux g, of the form: 



where D writes 



H 



dh 
dx 



D = 



(2.6) 



(2.7) 



and where the constant r (scaling with a volume) and the function / depend on the 
geometry and the flow equation and r] is the dynamic viscosity. Using the expression ( 2.6 1 



for the flux, eq. (2.3) admits a self-similar solution h{C,) 
variable C = x/^/Dt^ which obeys: 



H ipi^Q with the similarity 



This equation may alternatively be rewritten, in terms of C(^): 



2/ 



dip ^ dip 



(2.8) 



(2.9) 



The solution of the above equations can be found analytically or numerically, depending 
on the complexity of the normalized flux function f(ip). In the following, we will first 
recall the case of porous media, treated by |Huppert fc Woods(1995) and the 2D Stokes 
flow, addressed by Hinch(2007) (unpublished) and Taghavi et o/.(2009)Taghavi, Seon, 
[Martinez fc Frigaard[ We note that the latter paper included the effects of the rheological 
properties of the fluids. However, in order to focus on the geometrical aspects, we will 
assume in the following that both fluids are Newtonian and have the same viscosity. 

2.1. Lock- exchange in porous media 



For a homogeneous layer of porous medium of permeability k (see for instance Huppert 
I&: Woods(1995) ), the flow in each fluid is given by Darcy's law which relates the velocity 
in each phase to the local pressure gradient : 

KdP± 
rj dx 



(2.10) 



At a given location x, the velocity is then uniform in each layer, and the no net flux 
condition (eq. (2.5 )) simply writes: hux+ + {H—h)ux- = 0. The latter equation, combined 



with eq. (2.10) and cq. (2.2), leads to eq. (2.6), and thus (combined with eq. (2.3)) to eq. 
(2.8) with a diffusion coefficient and a flux function: 



DpM = K,H 



The solution of eq. ( 2.8 ), in the similarity variable C 



(2.11) 
(2.12) 

x/^DpMt, is then a linear profile 



(Huppert fc Woods(1995)): 



iP = h{x,t)/H= {1 + 0/2 



(2.13) 



The so-obtained front profile in homogeneous porous media is displayed in Fig. [2] (st raight 
line) together with the ones for rectangular cells (referred to in subsection |2.3[ ). The 
leading (i/i = 0, ^ = — 1) and traihng [ip — Q = 1) edges spread as y/DpMt. Therefore, 
the lock-exchange diffusion coefficient for porous media is DpM- It should be noticed 



that Bear(1988) reported a numerical integration of eq. (2.8) indicating that the gravity 
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Figure 2. Lock-exchange interface V'(C) between the two fluids for rectangular ceUs of different 
aspect rat ios F — H/b. The straight hne of slope 0.5 corresponds to the Darcy porous media 
limit (eq. ( |2.13| ). The other curves correspond to the case of a rectangular cross-section channel 
(subsection]^. From the bottom left to right, F = 10, 4, 2, 1, 0.5 and 0.2. 



current spreads as the square root of time. Note that the corresponding result for a Hele- 
Shaw cell, that is two parallel plates of height, H, separated by a tiny gap b [h <^ H), is 
obtained using the permeability k = b^/12: 



D 



HS 



b^HApg 
I2f] 



(2.14) 



2.2. Lock- exchange for a 2D Stokes flow between two horizontal boundaries 

For a 2D Stokes flow between two horizontal parallel boundaries, separated by a height 
H in the plane (z — x) (assuming invariance along the y— direction), the flow in each fluid 
is given by the Stokes equation: 



^2 / N 9P± 
rj V Ux±[x,z) = 



(2.15) 



At a given location x, the velocity profile consists of two parabola profiles matching 
the no slip boundary conditions at the bottom and the top boundaries (u2;+(a;, 0) = 
Ux-ix, H) = 0) and the continuity of the velocity {ux-{x,h) = Ux+{x, h)) and of the 
shear stress {rj dux-{x, h)/dz = rj dux+{x, h)/dz) at the interface, z = h. 



Using the no net flux condition (eq. (2.5|) and eq. (2.2), we obtain eq. (2.6), which enables 



to rewrite eq. (2.3), in the form of eq. (2.8) or eq. (2.9), with: 
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D 



H^Apg 
Si] 



(2.16) 
(2.17) 



Note that the polynomial development of the solution of eq. (2.9 1 around ip — gives: 
C = —Co + 2V'^/(3Co)- Thus, the location, (^{tp = 0) = —Co, of the leading edge of the 
interface is indeed constant in the similarity variable. Moreover, the development shows 
that the slope of the interface is vertical at the bottom wall {tp — 0). This is also the case at 
the upper wall {ip — 1), as the problem is symmetric with respect to the centre of the cell. 
We note that in the presence of such a vertical slope, our (horizontal) quasi-parallel flow 
assumption falls locally, but it is still valid upstream and downstream, where the slope 
of the interface remains small. The solution of eq. (|2.9[) can be found numerically using a 



shooting method similar to the one used by Hinch(2007) and Taghavi et a/.(2009)Taghavi, 



|Seon, Martinez fc Frigaard It was computed starting the integration of eq. (2.9) from 
{ip = 0.5, C = 0) and matching the asymptotic development in the vicinity oitp = 0. From 
the so-obtained solution, one can deduce the spreading diffusion coefficient between the 
leading edge {h ~ 0, —Co = —0.1607) and the trailing edge {h = H, Co) of the front, from 
[0.5 {x{h = 0) - x{h = H))]^ = Dt, which gives D2D = DQ, so that: 



D2D = 0.0086 i?' 



(2.18) 



This result is in agreement with the one found by |Hinch(20(37)] [Taghavi et aL(2009)Tagha\| i, 
Seon, Martinez & Frigaard provide five i/)(C) plots in their Fig. 9, corresponding to dif- 
ferent viscosity ratios and including our case. From that figure we may obtain a value 
of their similarity variable, 770 ^ 0.09, which is consistent with our finding Co = 0.1607, 
when taking into account their definition of the similarity variable, C = V VS- 

For completeness, D2D may be compared to the result for a cylindrical tube of diameter 



d (Seon et a/.(2007)Seon, Znaien, Salin, Hulin, Hinch & Perrin): 



Dt = 0.0054 



-3 Apg 

V 



(2.19) 



The above expression is indeed very close to the 2D result, with d playing the role of H. 

2.3. Lock-exchange in a rectangular cross-section channel 

This article aims to extend the computation of the lock-exchange diffusion coefficient to 
rectangular cells of arbitrary cross-sections H xb (see Fig.[T]). In the following, the cross- 
section aspect ratio will be denoted F = H/b. As in the previous section, a quasi-parallel 



flow approximation is assumed {i.e. small interface slope) which leads to eq. (2.2) for the 
pressure gradient. We will also assume the invariance of the interface location along the 
gap direction y. This requires that the deformation of the interface, induced by the flow 
profile along the direction y, relaxes much more quickly than the gradient along x. The 
flow in each fluid obeys a 3D Stokes equation: 



T] V^Ux±{x,y,z) 



dP± 
dx 



(2.20) 



In order to solve this equation, we follow the series decomposition in Fourier modes of 



the velocity field used by Gondret et a/.(1997)Gondret, Rakotomalala, Rabaud, Salin 



I&: Watzkyl This paper addressed the issue of the parallel flow of two fluids of different 
viscosities in a rectangular cell. This issue is very closed to ours, as it requires to solve the 
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Poisson equation (eq. (2.201), but with different viscosities and the same pressure gradient 



for both fluids in Gondret et aL(1997)Gondret, Rakotomalala, Rabaud, Salin & Watzky 



Poiseuille-like 



The method used was to split the velocity into two terms, Ux±{x,y,z) 

ul%{x,y,z). Here, the first term, u*±(a;,j/) = - (f )'] is the 

unperturbed velocity far away from the interface. The second term satisfies the Laplace 
equation, V^it*5-(x, y, z) = and vanishes far away from the interface. Its expression in 
terms of a sum of Fourier modes leads to a velocity profile of the form: 



^ 32(- 



62 dp+ 



Srj dx 
1)" (a±„e(2" 



1 



2y^2 



5±« e 



n=l 



(2n- l)^7r3 



■ cos[(2n — 1) 



'^y^ 



(2.21) 



in which the no slip boundary conditions at the two vertical walls {ux± {x^ y = ±6/2, z) = 
0) have been taken into account. Each Fourier mode, ((2n— l)7r/6), involves two constants 
for each fluid, a±n and b±n- These four constants are determined by using the no slip 
boundary conditions at the bottom and top of the cell {ux+{x, y,z = 0) = u^- (x, y,z = 
H) = 0) and the continuity of the velocity {ux+{x,y, z = h) = Ux-{x,y,z = h)) and 
of the shear stress (rj dux+(x,y, z = h)/dz = rj dux-{x,y, z = h)/dz) at the interface. 
Combining the so-obtained expressions for the velocity with eq. (|2.2|) and the no net flux 



(2.22) 



condition (eq. (2.5)), one obtains the horizontal flux of the heavy fluid (eq. (2.6)): 

dh 



-DHshih/H) 



with 



Dhs = 



/r(V') 



b^HApg 
1277 

-0 + ar(V') 



1 



7r 



(1 - V' - ar(V') - 7r) - Sr{ip) 



where 



ar(V') 



6rW 



7r 



96 (1 + e^(i"''')(2"-i)^) (1 _ e^''^(2«-i)7r-) 



r ^ 

n=l 



(l + er(2n-i)^) (2n- l)^7r5 



48(1 



„r(l-V) (2n-l)7r\ 



(1 



,r 4>(2n-l)Tr\ 



r ^ 

n=l 

CO 

n=l 



(-l + e2r(2n-l)7r) (2, 

192 tanh( ^('"~'^" ) 



l)^7r5 



(2n 



1)^5 



(2.23) 
(2.24) 

(2.25) 
(2.26) 
(2.27) 



Eq. (2.22) admits a self-similar solu tion, h{ () — Hip(C,), with the similarity variable 
C, = x/^Dhs t, which obeys eq. (2.8) or eq. (2.9). As previously, it is easier to compute 
the solution (^(tp) of eq. (2.9) subject to the corresponding asymptotics, C = —Co + 
8 ■03/(3^0) in the vicinity of the boundary, ip = 0. We solve this equation using the 
shooting method previously described and using Mathematica Software. The solutions 
h{() are plotted in Fig. [2]for different values of the cell aspect ratio F. We notice that, in 
contrast with Darcy predictions (straight line in Fig. [2]), but similarly to the case of the 
2D Stokes flow, the profiles, h{(), exhibit vertical slopes at the edges of the cell. We note 



also that such vertical slopes were observed in the experiments by Seon et aL(2007)Seon, 
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100 



r=H/b 



Figure 3. Log-log plot of the normalized lock-exchange diffusion coefficient, Dr/Dhs, versus 
the aspect ratio F = H/b, fo r rect angular cells. Dr is normalized by the diffusion coefHcient 
Dhs = HApg/ {12 Tj) (eq. (2.f4l, obtained for the Hele-Shaw cell limit, also corresponding 
to a 2D porous medium. Accordingly, the latter regime corresponds to the dashed straight line 
(Dr /Dhs = 1) in this representation. The solid straight line corresponds to the 2D Stokes limit 
(eq. (2.18l), leading to a slope 2 in this representation. 



Znaien, Salin, Hulin, Hinch & Perrin and Huppert & Woods(1995) When comparing 
their experiments in a Hele-Shaw cell with Darcy predictions, the latter authors reported 
that "Some discrepancies develop near the leading edge of the current as a result of 
the increasing importance of the bottom friction at the nose" (Fig. 2 of Huppert &: 
|Woods(1995)| ). This mismatch will be addressed in the next section. According to the 
so-obtained profiles, stationary in the similarity variable C, the leading and trailing edges 
of the front spread as the square root of time, and a lock-exchange diffusion coefficient, 
dependent on the cell aspect ratio, can be defined: 



b^HApg ,H, ,H, 



(2.28) 



Fig. [3] displays a log-log plot of the normalized rectangular cell lock-exchange diffusion 
coefficient, Dj^/Dhs = P{H/b), versus the aspect ratio F = H/b. At small aspect ratios, 
F < 1, the diffusion coefficient falls on top of the full line of slope 2, which corresponds to 
the 2D Stokes flow between boundaries distant of H {b 
ratios, F 



100, the diffusion coefficient approaches the dashed hue, Dr/Dhs 



00, eq. (2.18)). At large aspect 

1, 



obtained for the 2D homogeneous porous medium case (eq. (2.14)). We note that the 



latter case, which corresponds to a Hele-Shaw cell of infinite aspect ratio, overestimates 
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the lock-exchange diffusion coefficient, by a relative amount of about 30%, for aspect 
ratios as large as F = 10 — 20. 



3. 2D Stokes-Darcy model for lock-exchange in a rectangular 
cross-section channel 

The above-mentioned failures of the 2D Darcy model at finite aspect ratios may come 
from the velocity slip condition at the bottom and top edges of the cell {z = and 
z = H, respectively). This non physical condition is indeed required by the use of 
Darcy equation for the flow, which neglects the momentum diffusion in the presence 
of velocity gradients, in the plane of the cell {z — x plane). The momentum diffusion 
may however be taken into account in 2D, through the so-called Stokes-Darcy equa- 



tion (see Bizon et aL(1997)Bizon, Werne, Predtechensky, Julien, McCormick, Swift Sz 



Swinneyl Ruyer-Quil(2001) Martin et a/.(2002a)Martin, Rakotomalala & Salin Zeng 
et a?.(2003)Zeng, Yortsos fc: Salin[), which is similar to the Darcy-Brinkman equation 
used in porous media (see Brinkman(1947) ). This 2D model enables to handle discontinu- 



ities such as cell edges, gap heterogeneities and fluid interfaces (Ruyer-Quil(2001) Martin 



et aZ.(2002a)Martin, Rakotomalala & Salin Zeng et a/. (2003) Zeng, Yortsos & Salin Talon 



et aL(2003)Talon, Martin, Rakotomalala, Salin & Yortsos) and was successfully applied 



in the study of Rayleigh- Taylor instability (Martin et aL(2002a)Martin, Rakotoma- 



lala fc Salin[ Fernandez et aL (2002)Fernandez, Kurowski, Petitjeans & Meiburg Graf 



et aZ.(2002)Graf, Meiburg & Hartel), of dispersion in heterogeneous fractures (Talon 
et aL(2003)Talon, Martin, Rakotomalala, Salin & Yortsos I and of chemical reaction fronts 
(Martin et aL(20026)Martin, Rakotomalala, Salin & Bockmann). Although our present 
case of interest can be handled with 3D Stokes calculations, it is of interest to test the 
applicability of the 2D Stokes-Darcy model to the case of deep and narrow cells. Indeed, 
such a 2D model, once validated, could be a useful tool to address the issue of more 
complicated cases, such as gravity currents in the presence of viscosity contrasts, or in 
fractures with aperture heterogeneities. 

In this model, the flow in the rectangular cell (Fig. [T]) is assumed to be parallel to the 
plates {u{x,y,z) = {ux{x,y, z),0,Uz{x,y, z)) with a Poiseuille parabolic proflle across 
the gap (the key assumption). Using the Stokes equation with this y dependency, the 
gap-averaged fluid velocity U{x, z) — ^ J-b/2 ^(^' follows a Stokes-Darcy (SD) 

equation which reads here for the horizontal component of the velocity: 



^U.Mx,z) + /3 f^W^U^±{x,z) = ^ 



(3.1) 



to the Darcy's law (eq. 
mentioned above (eq. \2 



The flrst term on the le ft han d side of eq. (3.1) and the pressure gradient correspond 

j2 for the Hele-Shaw cell as 



2.10)) with a permeability k 



The second term on the left hand side of eq. (3.1) is 



the Brinkman correction to the Darcy equation (see Brinkman(1947) ), which involves 
an effective viscosity, /3r]. This effective viscosity may be taken equal to the one of the 
fluid (/3 = 1) for the sake of simplicity (or to enable the matching with a 2D Stokes 
regime at F 0). However, Zeng et aL (2003) Zeng, Yortsos & Salin showed that in 



the Hele-Shaw cell regime (at large F), the effective viscosity was slightly higher, with 
P = 12/7r^ ~ 1.215. At a given location x, integrating eq. (3.1) leads to the two velocity 



profiles matching the no slip boundary conditions at the bottom and the top boundaries 
{ux+{x,0) = Ux-{x,H) = 0) and the continuity of the velocity {ux-(x,h) = Ux+{x,h)) 
and of the shear stress {(3 ij dux-{x, h)/dz — (3r]dux+{x,h)/dz) at the interface. Us- 
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l.E-01 l£-hOO 1J;+01 1^+02 

r=H/b 



Figure 4. Log-log plot of the normalized lock-exchange diffusion coefficient, versus the aspect 
ratio r = H/b, for rectangular cells. The data points, obtained using the 2D Stokes-Darcy model 
(eq. ([4])) with Brinkman viscosity factors 13 = 1 (open squares) and /? = 12 /n'^ (open triangles) 
are compared to the full 3-D results already shown in Fig. [3] and depicted here as a solid line. 
The inset gives the relative difference between the 2D Stokes-Darcy model (/3 = 1, open squares 
and /3 = 12/7r^, open triangles) and the 3-D calculations. Note that these data were obtained by 
the difference between values of accuracy of the order of a few 10~^, which results in the small 
dispersion observed in the figure. 



ing the no net flux condition, Ux+dz + Ux~dz = 0, and eq. (2.2 1, we obtain the 
horizontal flux (eq. (2.4|): 



q{x) = ~DHs!sD{h/H) 



dh 
dx 



(3.2) 



where Dhs was already given in eq. (2.14) and the reduced flux function is equal to: 



fsu W = .... ' (2 + 4 (1 - ^) ^ - ^) + - 2 ^« 



where 



4d(rf-tanh(i) [ sinh(2d) 

(1 - ^) cosh(2 d (1 - ^)) cosh(2 d j;) _ cosh(d (1 - 2 V>)) ] 
+ sinh(2d) cosh(d) / 



and K — 6^/12. A comparison of the full 3-D calculations for a rectangular channel of 
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aspect ratio T = H/b with this 2D approximation can be performed on the flux functions, 



fri^) (eq- (2.24|) and fsnii') (eq- (3.3|). These two flux functions are close to each other, 
within a few per cents. In order to address the comparison in the range of interest for 
the Hele-Shaw assumption, i.e. F ^ 1, let us analyze the limit F — !■ cx) (d oo), which 
gives 

'3 ,,,, ,,,,./^ b , ^(,b.^ 



fsD,T^oo ^ ^{l - V') - (| - ^(1 - V'))y I + O [ij^r] (3.5) 
for the Stokes-Darcy flux and 

/r,r.oo ^(1 - ^) - (J - ^(1 - ^))l^^lf(^ 1 + O ((1)^) (3.6) 

for the full 3D rectangular cell flux (with Zeta{5) — J^T """^ ~ 1.03693, the value of the 
Riemann-Zeta function). We note that the leading term of both series corresponds to the 



expected porous media Darcy hmit (eq. (2.12)) with a permeability k = &^/12. However, 



the next order term {0{b/H)) is not the same, unless one chooses for the factor /3, 

/3 = 3 (^^^^^)2 ^ 1.192 (3.7) 



which is very close to the value 12/7r-^ ~ 1.215 found by Zeng et aL(2003)Zeng, Yortsos 



& Salin and to the value 6/5 proposed by Ruyer-Quil(2001) The lock-exchange diffusion 



coefficient has been computed, with the same procedure as above, by integrating eq. (2.9), 
using fsoitp), from ip — 0.5, and matching the asymptotics, ( — —Co + 8 ■0^/(9 Co) in 
the vicinity of the boundary, "0 = 0. The so-obtained lock-exchange diffusion coefficients, 
calculated for two different values of /3 (/? = 1 and f3 — 12/7r^) are compared to the 3D 
calculations in Fig. |4] The data for both values of /3 are indeed very close to the 3D 
data over the whole range of aspect ratios, F = 1 — 100. The inset of Fig. |4] gives the 
percentage of error for the two values of /3. We note that these data were obtained by 
the difference between values of accuracy of the order of a few 10~^, which results in the 
small dispersion observed in the inset of Fig. |4j As expected, the results at large F (in the 
Hele-Shaw regime) are closer to the 3D full problem for f3 = 12/7r^ than for /3 = 1. We 
point out that, whereas the Brinkman term does bring a significant correction, the exact 
value of the Brinkman viscosity factor /3 is however not crucial; For instance, for F = 10, 
we obtained a diffusion coefficient 3.5% smaller than the 3D value for /3 = 1 and 0.1% 
larger for /3 = 12/7r^, to be compared to the 30% of error if the cell was assumed to be 
of inflnite aspect ratio (Hele-shaw limit) as in Huppert fc Woods(1995)| In conclusion of 



this comparison, we have shown that the 2D Stokes-Darcy model for lock-exchange in a 
rectangular cell captures quite accurately the effect of the finiteness of the cross-section 
aspect ratio. By using the correct /? value, the error in the model is smaller than 5% for 
aspect ratios larger than F > 1. 



4. Experiments 

In this section, we will present experimental measurements of the diffusion coefficient 
in Hele-Shaw cells of different aspect ratios and we will compare them with our computed 
values. 

We used borosilicate rectangular cells of height H and thickness b and typical length 
30 cm (Fig. [l]). The rectangular cross-sections of the cells were (in mm^): 2 x 6, 2 x 12, 
2 X 20, 3 X 3, 3 X 9, 3 X 30, 4 X 6, 4 X 10, 6 X 6. Each cell was used with one side or the 
other held vertically, leading to two aspect ratios per cell. With such values, we covered a 
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Figure 5. Side view (plane z — x) oi the lock-exchange interface for rectangular cells of aspect 
ratios F = H/h = 10, 6, 2.5, 1.5, 1, 2/3, 0.4, 1/3, from top to bottom, respectively. The horizontal 
axis is scaled with V Dhs t, which allows to observe the decrease of i^o as F decreases. The vertical 
dashed line corresponds to C = 0. 

wide range of aspect ratios, from F = H/b = 1/10 to 10. We used, as Newtonian miscible 
fluids, aqueous solutions of natrosol and calcium chlorite. The fluids had equal viscosities, 
which were fixed by the polymer concentration and measured with an accuracy of 1%. 
The fluid densities were adjusted by addition of salt and measured with an accuracy of 
0.01%. The overall accuracy in D^s was typically 5%, when taking into account the 
above accuracies in viscosities and densities and the inherent temperature variations 
during the experiments. The viscosities and the densities of the fluids were chosen to 
satisfy two experimental requirements. The experiments must be fast enough in order to 
prevent any significant molecular mixing of the fiuids and one should be able to put the 
two fluids in contact without mixing. The latter condition requires a rather large density 
contrast and large viscosities. With our cell sizes, a good compromise was obtained with 
a density contrast of about 1% and typical viscosities in the range, 10 — 50 mPa.s, leading 
to a lock-exchange diffusion coefficient ranging from lO^^cm^ /s to \cw? js. The typical 
Reynolds number, built with the gap of the cell of these experiments is smaller than 0.1. 
For each experiment, the cell was first held with its axis Ox vertical. The fluids were 
successively slowly injected, with the lighter fluid on top of the heavier. Then the cell 
was closed and put in the desired position, with its axis Oz vertical in a few seconds. 
The development of the lock-exchange pseudo-interface was then recorded thanks to a 
video camera. Typical pictures (side view in the plane z — x) are given in Fig. [5] for 
cells of different aspect ratios. The horizontal axis is scaled with \/Dhs t, so that one 
can see the decrease of Co as F decreases. With this representation using the self-similar 
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Figure 6 . Plot of the square of the spreading distance of the leading edge of the front versus 
time for five cells: (•) H = 3 mm, b = 30 mm; (□) H = 6 mm, b = 2 mm; (*) H — 6 mm, 
b = 4 mm; (o) H = b = 6 mm; (+) H = 30 mm, b — 3 mm. The solid line is a linear fit to the 
data, the slope of which gives the diffusion coefficient plotted in Fig. IT] (top). 



variable C, — x/ \/Dhs the profiles are stationary. One may notice that the trailing 
edge is fuzzy. This can be attributed to the stick condition at the upper wall: The dark 
dense fluid does stay at the walls for a long time, in particular in the corners of the cross- 
section. The same phenomenon takes place at the bottom of the cell, but the presence 
of transparent light fluid has little effect on the turbidity of the heavy dark fluid, and is 
therefore not noticeable on the pictures. It is worth noting that the shape of the leading 
edge evolves from an edge at large aspect ratios F to a more and more step-like shape 
as F decreases. It should be noticed that for small aspect ratios, although it is rather 
difflcult to take pictures, a top view of the cell reveals a mild spanwise dependency 
of the interface, but we do not observe the spanwise lobe-and-cleft instability reported 
by Simpson(1972) For each experiment, the locations of the leading and trailing edges 
of the front were measured in time. Fig. [6] gives the variations of the square of the 
spreading distance versus time for five cells. It is worth noting that the dependency 
is almost linear: Therefore a linear fit provided the lock-exchange diffusion coefficient, 
with a typical accuracy of 20%. Fig. [t] (top) displays the so-obtained normalized lock- 
exchange diffusion coefficient as a function of F. One can see that the agreement with the 
3Z? calculations over the two decades of our measurements is rather good. We note that 
for the large aspect ratio limit of the experiments (up to F = 10), the Hele-Shaw cell limit 
is not reached, and would underestimate, by 30%, the lock-exchange diffusion coefficient. 
This result thus confirms that for such aspect ratios, one should either compute the 
full 3D Stokes equation or use the Stokes-Darcy model to obtain the correct behaviour. 
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Figure 7. Top: Log-log plot of the normalized measured lock-exchange difltusion coefficient, 
Dfixp/ Dhs (open squares), versus the aspect ratio, F = H/b, of the cell. The solid line cor- 
responds to the 3-D results (same as in Fig. |4|. Bottom : Superimposition of the measured 
pseudo-interface between the two fluids (grey fuzzy line) and of the theoretical profile (full dark 
line) corresponding to Fig. [2] The aspect ratio is F = 4. 



We also note that our calculation still holds for aspect ratios as small as F = 0.1. This 
result is quite unexpected since for such aspect ratios, some spanwise dependency of the 
profile was observed, and the hypothesis of the interface surface, h{x,y,z), invariant in 
the y direction is certainly broken. The bottom of Fig. [7] displays the superimposition 
of the theoretical and the experimental interfaces between the fluids, for an aspect ratio 
F = 4. The agreement between the two is rather good, thus validating our model. Such 
an agreement is rather surprising as our small slope assumption is violated at the edges 



of the gravity current. This agreement, already emphasized by Huppert(1982) and Seon 



et aL (2007) Seon, Znaien, Salin, Hulin, Hinch & Perrin is likely to be common to viscosity 
dominated gravity current without surface tension. 



5. Conclusion 

The viscous lock-exchange diffusion coefficient refiects the competition between the 
buoyancy driving effect and the viscous damping, and depends on the geometry of the 
channel. We give the backbone to calculate this coefficient in different configurations: We 
recall its computation for a porous medium already found by |Huppert fc Woods(19957 



and compute it for a 2D Stokes flow between two parallel horizontal boundaries separated 
by a vertical height, H. This result is in agreement with Hinch(2007)|(unpublished) and in 
reasonable agreement with recent computations by Taghavi et aL(2009)Taghavi, Seon, 
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[Martinez fc Frigaard l Using a quasi-parallel flow assumption, we have calculated the 
pseudo-interface profile between the two fluids and the diffusion coefficient of viscous 
lock-exchange gravity currents for a rectangular channel (horizontal thickness 6, vertical 
height H) of any aspect ratio {H /b). This analysis provides a cross-over between the 2D 
Stokes flow between two parallel horizontal boundaries separated by a vertical height, H, 
and the Hele-Shaw cell limit (applying for H/h > 100). Moreover, the shape of our profiles 
allows to account for the discrepancy observed at the nose of the gravity current in the 
experiments by Huppert fc Woods(1995), The agreement, obtained despite the failure of 
the lubrication assumption at the edges of the current, should deserve however further 
theoretical investigation. Our calculations of the diffusion coefficient and of the shape of 
the profile have also been convincingly compared to new experiments carried out in cells 
of various aspect ratios (1/10 — 10). We have also calculated the lock-exchange diffusion 
coefficient for the same rectangular cells, using the 2D Stokes-Darcy model. This model 
is shown to apply to aspect ratios H/b > 1, provided that the appropriate Brinkman 
correction is used. Such a 2D model may be useful to describe gravity currents with a 
finite volume of release, with fluids of different viscosities, or in heterogeneous vertical 
fractures. 
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